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Realizing a quantum transistor built of molecules or quantum dots has been one of the most 
ambitious challenges in nanotechnology. Even though remarkable progress has been made, being 
able to gate and control nanometer scale objects, as well to interconnect them to achieve scalability 
remains extremely difficult. Most experiments concern a single quantum dot or molecule, and they 
are made at ultra low temperature to avoid decoherence and tunneling. We propose to use canonical 
transformations to design quantum devices that are protected by symmetry, and therefore, may be 
operational at high temperatures. We illustrate the idea with examples of quantum transistor 
architectures that can be connected both in series and parallel. 

PACS numbers: 85.65.+h, 73.63.-b 



I. INTRODUCTION 

Thanks to advances in nanofabrication, experimental- 
ist can routinely manufacture nanostructures that resem- 
ble artificial atoms -quantum dots- that can be manipu- 
lated with an extreme degree of control^. The degree of 
tunability is such, that one can see a state with a single 
electron sitting at the dot, and even control and detect 
its spin. By tuning a gate voltage one can also control 
the transport through the dot in what is regarded as a 
single-electron transistor. The techniques used to con- 
struct these devices can also be applied to build struc- 
tures with multiple dots, and suggest that scalability 
to circuits of many quantum dots should be technolog- 
ically achievable. Another avenue that has been pur- 
sued is using molecules to build transistors, instead of 
quantum dotsi Some examples of working devices have 
been realized with nanotubes^, benzene molecules^—, 
fulleren o 10 ! 11 , graphen o 12 ' 13 , and proteins^, to mention 
a few. Adding contacts and gating these devices is very 
challenging, and scalability may be more difficult to at- 
tain. A third possibility consists of manipulating single 
atoms to build complex structures^ that could be con- 
trolled using an STM tipi£. 

Whether a complex molecule, or a quantum-dot ar- 
ray, the most general transport features of these sys- 
tems in the absence of interactions can infered from 
some basic fundamental physics: Aharonov-Bohm-typc 
interference^, and Fano resonances^. There is an 
important effort toward engineering devices at the 
atomic/molecular scale taking advantage of these differ- 
ent phenomena, and that can have useful functionality 
as classical switches -diod, transistor-, or novel behav- 
ior that could open the doors to quantum circuits with 
no classical analogy. These systems could lead to a new 
generation of electronic devices, and revolutionary appli- 
cations such as quantum computers. 

The chemistry of these molecules, or the geometry of 
the circuits, can be simple, or not. In this work we show 
that regardless of the chemistry, the transport proper- 
ties of a molecular device can be inferred from symmetry 
considerations^^. Moreover, canonical transformations 
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FIG. 1: Illustration showing the original system (a), and the 
equivalent one after the folding transformation is applied to 
the sites 0, ...L — 1 (b) and (c). The emerging n phase can be 
eliminated by starting from the circuit (d). The value of the 
hopping connecting the two halves becomes t* = t/\ /f 2 after 
the transformation. 



allow us to establish dualities between quantum circuits, 
providing us with a tool to design devices that are pro- 
tected by symmetry. 



II. EQUIVALENT QUANTUM CIRCUITS 

Consider a one dimensional chain of atoms, described 
by a tight-binding chain with hopping t. with a weak link 
in the center, with hopping t q . Let us split the chain into 
left and right halves, and fold it in two, as depicted in 
Fig. [Ha). The corresponding Hamiltonian is: 
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where c\ t j is the electron annihilation operator acting 
on site j of lead A (where the values 0,1 correspond to 
left and right leads), and sites labeled "0" connect the 
two leads. The addition of the term proportional to V g 
will become apparent below. For the moment, we take 
V g = 0, and consider spinlcss fermions, but these consid- 
erations can be easily generalized to the spin-full case. 
We now introduce a symmetric (+) and antisymmetric 
(-) combination of operators acting on the left and right 
leads. This is nothing else but an application of the re- 
flection symmetry, yielding new even(+) and odd(-) op- 
erators: 
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This is a simple change of basis, with the new operators 
obeying the same fermionic anti-commutation rules. This 
is a well known canonical transformation, originally in- 
troduced in the context of quantum impurity problems^, 
and referred to as a "folding" transformation^. It is used 
extensively in Numerical Renormalization Group (NRG) 
calculations^, and it has been recently used to study im- 
purity problems on rings using the density matrix renor- 
malization group (DMRG)2i. It has also been used in 
rather creative ways to simplify complex problems in 
quantum chemistry, such as conjugated molecules and 
polymers^. In this work, we shall take an unconven- 
tional approach: instead of applying the transformation 
to the entire system, we shall apply it to only a fraction 
of it, say L sites to the right of the weak link. As a result, 
we find a section of the chain transformed to the new ba- 
sis, connected to the rest of the left and right leads that 
remain in the original basis, as shown in Fig. 1(b). At 
the boundary connecting the transformed and untrans- 
formed sections, we find new hopping terms mixing the 
two bases as: 
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After some simple algebra, it is easy to see that the hop- 
ping at the weak link, maps onto effective chemical po- 
tentials in the new basis (see Appendices A and B for 
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FIG. 2: Illustration showing the the transformation that 
leads to an annulene-like geometry: (a) original system, and 
the equivalent one after the folding transformation is applied 
to the sites 0, ...L — 1 (b) and (c). The emerging n phase can 
be eliminated by starting from the circuit (d). The arrow in 
panel (a) indicates a dangling link. 
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where n± = c < j s .c± is the density operator. This term is 
just a boundary chemical potential with opposite signs 
for the "+" and "-" sections of the leads, at site "0" . All 
the other hopping terms along the "+" and "-" sections 
of the lead remain unaltered in the new basis. 

Taking all this into account, and after some visual in- 
spection, it is easy to see that the original chain, after 
the basis transformation, maps onto the equivalent cir- 
cuit shown in Fig. HJc): a diamond connected to the left 
and right leads (that remain in the original basis), and 
the "+" and "-" leads that have a finite length, and a 
boundary chemical potential at the ends of magnitude 
t q . Let us call these leads "gates" from now on. All 
hoppings inside the diamond are equal to t/y/2, with the 
exception of one link that has the opposite sign. This 
can be interpreted as a i magnetic flux threading the 
diamond. 

The physics of the original system, and the new equiv- 
alent one should be completely equivalent, since this 
mapping is exact. Therefore, let us consider setting the 
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III. CANCELING THE MAGNETIC FLUX 
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FIG. 3: Rearrangement of the atoms in the circuit: (a) single 
tight-binding chain and (b) T-junction. In both cases a por- 
tion of the system is hanging, leading to an interference effect 
due to a Fano resonance, as depicted in the inset. 



boundary chemical potential of the two "+" and "-" leads 
on the equivalent circuit to zero, t q = 0. If we look at 
the original chain, this corresponds to effectively cutting 
the connection between left and right leads, completely 
opening the circuit. Therefore, under these conditions, 
no current can circulate from left to right leads, no mat- 
ter what the temperature is. The circuit is completely 
open, by simple symmetry considerations! 

If one wanted to describe the transport properties of 
the equivalent circuit, one could use ideas of quantum 
interference^—, and would consider the anti-resonances 
coming from the gates^S or combinations of the two2i. 
However, this simple mapping shows a different aspect of 
the problem, and allows one to understand all the trans- 
port properties by simply studying at the original un- 
folded tight-binding chain. 

So far, we have shown that the equivalent circuit acts 
like a multi-gate transistor: by simply tuning the poten- 
tial at the gates, we can open the circuit, and completely 
suppress transport between the leads. Alternatively, at 
finite tunneling amplitude tq ^ 0, we can obtain per- 
fect transport. Another situation that we can consider is 
connecting one of the gate leads to ground, and the other 
one to a gate voltage V g . This corresponds to choosing 
t q = 0, V g ^ in Eq.®: 
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The results for the conductance are shown in Fig|4j At 
zero voltage V g = 0, the circuit is open, while at large 
voltages we recover the perfect conductance. This be- 
havior is analyzed below, and in Appendix C. 



One of the undesired aspects of the previous "device" is 
the emergence of a 7r flux threading the loop. One could 
easily devise a setup to get rid of the phase. Consider the 
circuit shown in FigJTJd), corresponding to a T-shape 
junction. After applying the folding transformation to 
all the sites situated to the left of the connecting link, we 
recover the same molecular circuit with the diamond in 
the center, but with no phase shift. The price we have 
to pay is that we can no longer interpret the problem in 
terms of a simple tight-binding chain. However, this still 
is a non-interacting circuit, and we can easily understand 
the transport properties. In this case, the conductance 
will depend on the parity of the dangling lead in the T- 
junction. Let us first consider the V g = case. If t q ^ 
this segment will have an even number of sites, with a 
vanishing density of states at u> = 0. This means that it 
cannot participate in the transport between the left and 
right leads, yielding a perfect conductance. If t q = 0, the 
segment is effectively cut in half, and now two situations 
can arise: (i) for L = Ak, we obtain the same result as 
before for t q ^ 0, or (ii) for L — Ak + 2, a resonance 
will appear at the Fermi level, and the conductance will 
cancel identically. Notice that these are more complex 
generalizations of the problem studied in Ref»22. The 
general situation for finite V g will be studied below. 



IV. GENERAL CASE: ANNULENES 

We can generalize the previous protocol to other ge- 
ometries, to obtain the equivalent of molecular transis- 
tors with a ring structure of N sites, similar to the pro- 
posed devices based on annulenes^. Consider the geom- 
etry illustrated in Fig 121 a). It depicts a one dimensional 
chain, that has been folded in half at the weak link, but 
with the addition of a second "twist" plus as dangling 
link, marked with an arrow in Figj2fa). By rotating the 
(N — l)/2 operators to the left of the Lth site, we obtain 
the hexagonal molecule shown in Figl^b), connected to 
leads, and gates, as in the original case described be- 
fore, with the emerging tt flux. The main difference with 
the tight-binding chain is the dangling sites, which in- 
troduce an additional degree of complexity. If the size 
of the annulus is N — Ak, the dangling link will have 
an even number of sites, and will not participate in the 
transport properties, yielding the same results as the di- 
amond. For N = Ak + 2, the dangling link will introduce 
an anti-resonance, and the conductance will be G = 0. 

Again, the undesired n phase can be removed by re- 
shaping the circuit as a T-junction (see Figj2jd)). In this 
case, we expect a parity effect in the transport behavior, 
depending on the size of the ring A, and the size of the 
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V. TRANSPORT PROPERTIES 

We now analyze the case with t q = 0, V g ^ in more 
detail. In general, the transport properties are related to 
the Green's functions that propagate one electron from 
the left to the right leads Gl_r(w). The conductance G 
at zero bias can be written as^i 

G=^t 4 \G LR (E F )\ 2 Po (E F ) 2 , (6) 

where p is the local density of states (LDOS) at the 
leads (assumed to be semi-infinite tight-binding chains) 
and Ep = is the Fermi level. 

Since the system is non-interacting, we can calculate 
the exact Green's functions using either the equations of 
motion formalism^ or by applying Wick's theorem. 

In Figure [3] we have rearranged the geometry of the 
leads to show more clearly the structure of the system. 
Independently of the phase, the system always presents 
a dangling segment. The conductance is controlled by 
Fano resonances. Two different paths, as represented in 
the inset of Fig. [3{a) , can interfere in a constructive or 
destructive way2&. The parity of the number of sites in 
the ring will determine whether the phases cancel or add 
up. 

It is convenient to start the discussion by considering 
the T-junction problem. In this case, the transport prop- 
erties depend on both, the parity of TV, and the size of the 
gate leads L g . We first look at a ring with N = Ak sites. 
In Figure 2] we present the results for the conductance 
for rings with 4 and 8 sites as a function of V g and gate 
leads with different lengths L g , odd and even. We find 
two clearly distinct situations depending on the parity of 
L g . As shown in this figure, when L g is even a perfect 
conductance of G = 2e 2 /h independent of V g is realized, 
while for odd L g the conductance depends strongly of 
the gate potential V g . To understand this behavior we 
look at Figure EJa). As explained before, in the V g = 
limit the upper part of the system is disconnected and 
the conductance is controlled by the number of sites in 
the gate leads L g . If L g is odd, this segment will have a 
resonance at the Fermi level and again, this will cancel 
the conductance identically. If L g is even, there will not 
be any resonance in the LDOS and the conductance in 
this case is perfect, G=2e 2 /h. 

When V g ^ the upper part of the system of Fig- 
ure [3ja) will be connected and the chemical potential at 
the ends will also start playing a role. The total num- 
ber of sites of the hanging part will be 2L g + (N — 2)/2 
always odd, producing, in principle, a resonance at the 
Fermi level, and canceling the conductance. However, at 
large V g the chemical potential removes one site from the 
hanging segment, with the conductance growing toward 
2e 2 /h with increasing V g . 

In order to understand the effects of a finite gate volt- 
age V g , let us consider the system before the transforma- 
tion, as shown in Figures [TJa) and (d), and take L = 1, 
such that we apply the folding transformation only on the 



link connecting left and right leads i?imk ©. This will 
yield two orbitals "+" and "-" that will be connected to 
the rest of the system through i/boundary ©■ In the large 
V g limit, the orbital "+" will be empty, and effectively re- 
moved from the system, which will have one site/orbital 
less. This changes the parity of the problem, converting 
a resonance into an anti-resonance and viceversa. 

We now analyze rings with N = Ak+2 sites. Results for 
the conductance as function of V g for rings with 6 and 10 
sites are shown in Fig. [SJ Same as before, we can observe 
two different regimes. For L g odd the conductance is 
always zero while for L g even the conductance depends 
strongly on V g . The physics can be explained following 
similar arguments as the previous case, but with odd and 
even L g cases inverted. 

We now go back to the case with a flux, corresponding 
to a single chain. As shown in Figures 0] and [SJ the con- 
ductance is independent of L gi as we should expect for 
a tight-binding chain. However, there is a parity effect 
corresponding to the size of the ring being N — Ak or 
N = Ak + 2, which dominates the transport properties. 
As mentioned before, this is explained by the effect of the 
dangling links added to the chain to construct the rings. 
Depending on the number of sites in the links, we will 
observe a resonance or anti-resonance at the Fermi level 
(see FigEPO). 

The effect of V g can again be explained as gradually 
removing a site from the chain (Appendix B). However, 
since adding or removing a site from a chain does not 
have any effect on its fundamental transport properties, 
the conductance interpolates smoothly between the t q = 
0, V g = and t q ^ 0, V g — cases. 

For both TV = Ak, and N = Ak + 2 cases, with L g odd, 
the behavior is exactly the same for the chain and the 
T-junction. This is a remarkable result, and leads one to 
think of a hidden symmetry between the two setups. As 
a matter of fact, it is not a symmetry relating the two 
Hamiltonians, but a symmetry in the Green's functions 
that determine the equilibrium transport properties, that 
is "magically" realized when these quantities are evalu- 
ated at the Fermi level, uj = 0. This is analyzed in great 
detail in the Appendix C. 



VI. MOLECULE IN PARALLEL 

Let us look again at FigQJd). We can make the trans- 
formation inside the hanging segment of the T-junction, 
away from the contact to the leads, as shown in FigJHJa). 
It is easy to see that the equivalent circuit results as de- 
picted in FiglHKb). The physics of the problem, of course, 
is just as described in the previous sections. It is the same 
system after all! This shows that connecting the molecule 
in parallel, or in series, will have exactly the same effect 
and functionality. Moreover, it may be easier to realize 
experimentally, and more convenient for scalability pur- 
poses. 
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FIG. 4: Conductance for rings with 4 and 8 sites (iV = 4fc 
sites), and rj 9 = 0. Two regimes appear: when L g is even, the 
conductance is perfect 2e 2 /h, while for L g odd, it strongly 
depends on V g . 
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FIG. 5: Conductance for rings with 6 and 10 sites (N = 4fe + 2 
sites) and t q = 0. Again, two different regimes are realized: 
When L g is odd the conductance is always zero. For L g even, 
the conductance depends strongly on V g . 



VII. SUMMARY AND CONCLUSIONS 

We have illustrated how three complex quantum "de- 
vices" , shown in Figures QJc) , f^c) , and EJb) can be un- 
derstood in terms of simple equivalent circuits that have 
the general topology of either a single chain, or a T- 
junction. These exact mappings have been obtained as 
a straightforward application of the folding transforma- 
tion. We can easily understand the transport properties 
of these simple circuits, and explain the general behav- 
ior of the molecular device. Two remarkable aspects are 
that: (i) in the single chain case, with a phase n threading 
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FIG. 6: Illustration showing the equivalence between (a) a 
T-junction, and (b) a molecular device connected in parallel. 
The equivalent circuit without a phase is depicted in (c). 



the molecule, the gate voltage can be used to effectively 
"cut" the equivalent circuit in two, yielding a complete 
cancellation of the conductance. This effect is protected 
by the symmetry of the problem, and indicates that this 
behavior as a quantum switch survives at arbitrary tem- 
peratures, (ii) The zero-temperature transport proper- 
ties of the chain and the T-junction are the same for L g 
odd, and can be understood in terms of a symmetry in 
the equilibrium Green's functions at lu = 0. 

Even though we have limited our analysis to a sin- 
gle orbital situation, these considerations can be gener- 
alized to a multi-orbital case. The non-interacting limit 
assumed here applies to a problem that can be described 
in terms of a Huckel-like theory. The addition of inter- 
actions would introduce long range terms in the trans- 
formed basis, making the folding transformation imprac- 
tical. However, one could imagine studying the inter- 
actions within a Hartree-Fock framework, in which the 
resulting effective Hamiltonian is quadratic. In this pic- 
ture, the interactions introduce a shift in the diagonal 
site energies, and the main effect is to shift and to split 
the resonance position o 37 ' 38 . In such scenario, as long 
as the reflection symmetry is preserved, our treatment is 
still valid. 

Inteference-based molecular transistors have already 
been propose d 27 ' 39 ' 40 , and quantum interference has re- 
cently been observed in molecules 4 ^. This study indicates 
a path toward understanding and designing molecular 
devices with customized functionality, and we illustrate 
it with setups that operate as multi-terminal transistors 
that can work connected in series, and parallel. These 
symmetries make us see these systems in a new light, 
and can also aid ab-initio calculations. It is clear that 
the mapping is not limited to the particular examples 
presented here, but it is more general and may lead to 
complex and interesting possibilities that could pave the 
way toward new advances in nano-electronics. 
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FIG. 7: Diagrams showing the resulting circuit after the 
transformation is applied to the encircled sites on the left. 



The interesting situation arises when the transforma- 
tion is applied to the left of certain boundary, and not to 
the sites to the right. Then, at the link connecting the 
transformed sites, and those remaining in the original 
basis, new hopping terms emerge. The four possibilities 
shown in Figure EJc)-(f) correspond to: 



(11) 
(12) 

(13) 
(14) 



By simply combining these diagrams, and introducing 
"twists" in the chains, we can construct the devices pre- 
sented in the body of this work. Of course, this is just a 
limited number of the countless possibilities. 
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APPENDIX A: SYMMETRY 
TRANSFORMATIONS 

The entire treatment utilized in this work relies on a 
simple and basic symmetry : reflection, equivalent to a 
bonding/ anti-bonding transformation. 
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c o/U = ±c_j). (7) 



We start by considering a non-interacting one- 
dimensional chain, folded in two. It is easy to under- 
stand the effects of this symmetry by visualizing them 
in diagrams, as depicted in Figure [7] When this trans- 
formation is applied to two parallel links, as in diagram 
[7ja), the transformation leaves the Hamiltonian invariant 
and reads exactly the same in the transformed basis: 
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When applied to a link connecting a site on chain "0" 
and its partner on chain "1" , it leads to a vanishing hop- 
ping, and the introduction of two chemical potentials (see 
diagram (7{b)): 
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(9) 
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In this Appendix we give a closer look at the link term 
connecting the two chains: 
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(15) 

After the considerations introduced in the previous Ap- 
pendix, it is easy to see that this Hamiltonian can be 
transformed into: 
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When t q = 0, this reduces to: 



-Hiink = V„n 



(16) 
(17) 



For large \V g \, this term will force the "+" orbital to be 
either fully occupied, or empty, with the consequent sup- 
pression of charge fluctuations. Therefore, all the trans- 
port will occur through the "-" orbital, meaning that the 
total system will behave as though it had one less site. 



APPENDIX C: SYMMETRIES IN THE 
CONDUCTANCE 

In order to understand why the conductance is the 
same for the T-junction and the single chain, i. e. does 
not depend on the flux for N — 4fc, and N = Ak + 2, 
L g odd, we need to calculate the propagator Glr ap- 
pearing in the expression for the conductance. We first 
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realize that, independently of the connection, both sys- 
tems have a common part shown in Fig.|SJa). Therefore, 
we calculate the partial propagators for this common seg- 
ment, g a R and gpR, and we establish the connection to 
the second lead using a Dyson equation. 

For the T-junction, the left lead must be connected to 
the site labeled as a, while for the the single chain, the 
connection must be established through site j3. This is 
shown in Figures H£b) and (c). 

The equations of motion for the propagator Glr for 
the T-junction case can be written as, 
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while for the chain we have, 
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where t^ a = t^p = t, and we had preserved the subindcx 
to distinguish the different cases. The bare Green func- 
tions gij are the propagators before reestablishing the 
connection between sites a and R. The Green's func- 
tions gij denote the propagators after reconnecting sites 
a and R as shown in Fig. Uta). 

Solving these systems of equations, we obtain for the 
first, 
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and from the second set of equations, 
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(21) 



(22) 
(23) 



Comparing Eqs. (|20)) and (|21j) . we can observe that 
the term between brackets is common to both cases. The 
prefactors will determine when these two Green's func- 
tions are equal. Note that by symmetry g aa = gpp- The 
prefactors in Eqs. (f20|) and (|2T|) as well as g aa and gpp 
from Eqs. ([22]) and ([23]) will be equal if, 



Qoccx — Q$ql- 



(24) 



Notice that, since we are calculating the zero bias con- 
ductance, this condition must be proved just for u> = 0. 

In order to demonstrate that the condition (f2"4")l is satis- 
fied, we need to calculate the bare Green's functions for a 
single tight-binding segment of length n = 2L g + (N — 2). 
For simplicity we will restrict this treatment to rings with 
N = 4 sites, n = 2L g + 2, but the demonstration can be 
extended to arbitrary ring sizes. The Green's function is 



FIG. 8: Rearrangement or the different components of the 
circuit: (a) the common block for both systems. Connections 
for the T-junction (b) and for the chain (a) systems. 



the inverse of the matrix (H—Iz) where H is the Hamilto- 
nian of the segment. For the non-interacting case studied 
here, H is a finite n x n symmetric tridiagonal matrix: 
Hij = cii for i — j, and Hij = bi for i = j ± 1. There is 
a simple and elegant formula to calculate the inverse of 
such matrix^: 

gij = (-l) i+j bi ■ ■ ■ bj-i9i-i<f) j+1 /9 n for i < j, 
gij = (-l) i+J 6j • ■ ■ bi-xOj-^i+x/On for i > j, 

where the O^s verify the recurrence relation 

6> 4 = diOi-i - 6?_ 1 ^i_ 2 , for i = 2, • • • , n, 

with 9q = 1 and 9\ = oi, and the </>i's obey the recurrence 

<f>i = a l (j) l+ i - b 2 (j)i + 2, for i = n - 1, • • • ,1, 

with (/> rl +i = 1 and 4> n = a n . Notice that 9 n = det(H — 
Iz). 

We are interested in the Green's function at the Fermi 
level, z = (we have ignored any imaginary contributions 
for simplicity). Moreover, all the diagonal elements are 
a; = except a n /2 = a n /2+i = V g /2. Similarly, all the 
off-diagonal elements are bi = —1, with the exception of 
Vg/2. Using these expressions, it is easy to verify 



»n/2 

that 9 n = dct(H) = for n = 4k + 2, and 
1 for n = 4fc. Moreover, for n = 4k, 



det(H) = 



g nl = (-l) 1+n b 1 ---b n ^ 1 =-V g /2, 



and 



fh 



_ ( i \n/2+li2 .2 . 

i - o n _ 2 b n _ A ■■■b 



l n/2°n/2°ri/2-2 ' 



b\b. 



2,2 



-Vg/2. 



Therefore, we can see that the recursions lead to the 
same numerical result for gn and g n \ as long as the 
length of the chain obeys n — 4k = 2L g + 2, or L g 
odd. The important condition to get this equivalence 
is the presence of the symmetric chemical potentials V g 
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in Hunk- Notice that these observations do not imply at the Fermi level and, consequently, their equilibrium 
a symmetry relating the Hamiltonians of the chain and transport properties, 
the T-junction, but a symmetry in the Green's functions 
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